Stochastic Neutrino Mixing Mechanism 
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We propose a mechanism which provides an explanation of the Gallium and antineutrino re- 
actor anomalies. Differently from original Pontecorvo's hypothesis, this mechanism is based on 
the phenomenological assumption in which the admixture of neutrino mass eigenstates in the mo- 
ments of neutrino creation and detection can assume different configurations around the admixture 
parametrized by the usual values of the mixing angles 0i2, O23 and 813. For simplicity, we assume a 
Gaussian distribution for the mixing angles in such a way that the average value of this distribution 
is given by the usual values of the mixing angles and the width of the Gaussian is denoted by a. 
We show that the proposed mechanism provides a possible explanation for very short-baseline neu- 
trino disappearance, necessary to accommodate Gallium and antineutrino reactor anomalies, which 
is not allowed in usual neutrino oscillations based on Pontecorvo's original hypotheses. We also 
can describe high-energy oscillation experiments, like LSND, Fermi and NuTeV, assuming a weakly 
energy dependent width parameter, a{E), that nicely fits all experimental results. 

PACS numbers: 14.60.Pq 
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I. INTRODUCTION 



Neutrino oscillations are a direct consequence of the 
assumption raised in the seminal article by Bruno Pon- 
tecorvo in 1957 jT] which asserts that neutrino states 
interacting with charged leptons through weak inter- 
actions are superpositions of neutrino states of non- 
vanishing definitive mass. In his paper, Pontecorvo used 
an analogy with neutral kaon mixing to propose that 
neutrino-antineutrino transitions may occur. Although 
such matter-antimatter neutrino oscillation has not been 
observed, this idea formed the conceptual foundation 
for the quantitative theory of neutrino flavor oscilla- 
tions, which were first developed by Maki, Nakagawa, 
and Sakata in 1962 [2] and further elaborated by several 
authors [MT]. 

Under the assumption that flavor eigenstates are dif- 
ferent superpositions of mass eigenstates, neutrino fla- 
vor oscillations can be described in the following way: 
as a neutrino propagates through space, the quantum 
mechanical phases of the mass eigenstates advance at 
slightly different rates due to the tiny differences in the 
neutrino mass eigenvalues. This results in a changing 
admixture of mass states as the neutrino travels. But a 
different admixture of mass states corresponds to a dif- 
ferent flavor state. So a neutrino born as, say, an elec- 
tron neutrino will be some different admixture of elec- 
tron, muon, and tau neutrino after traveling some dis- 
tance. Since the quantum mechanical phase advances in 
a periodic fashion, after some distance the state will re- 
turn to the original admixture, and the neutrino will be 
again an electron neutrino. The electron flavor content 
of the neutrino will then continue to oscillate as long as 



the quantum mechanical state maintains coherence. It is 
because the mass differences between the neutrinos are 
small that the coherence length for neutrino oscillation is 
so long, making this microscopic quantum effect observ- 
able over macroscopic distances. 

Interesting enough, neutrino flavor oscillations are the 
basis of the so-called solutions of several puzzling neu- 
trino observations recorded along the last four decades. 
The solar neutrino deficit initially observed in different 
experiments, counts on the neutrino oscillations reso- 
nantly enhanced by solar matter, the MSW phenom- 
ena [H in], to explain the solar neutrino observations. 
Similarly the same oscillation parameters can explain the 
deficit in Kamland experiment. Also the strong zenith 
dependence of atmospheric neutrino and antineutrino 
data can be explained by evoking neutrino oscillations. 
Finally, completely different sources of muon-neutrinos 
and muon-antineutrinos produced by meson decays in 
accelerators, confirms the necessity of neutrino oscilla- 
tions to understand the observations [TU|. Furthermore, 
recent measurements of experiments collecting neutrinos 
from reactors observed the necessity of a nonvanshing 
neutrino mixing ^13 [TI] [T^] , composing a robust picture 
in favor of the Pontecorvo's hypotheses which give rise 
to neutrino oscillations. This complete scenario involv- 
ing the three neutrino generations was recently analyzed 
intheRef. ^3\. 

On the other hand, one can argue that both Pon- 
tecorvo's hypotheses, i.e. neutrinos are massive and there 
exists neutrino mixing, have been experimentally con- 
firmed only indirectly through their main consequence, 
precisely, the neutrino quantum oscillations. In fact, the 
first Pontecorvo's hypothesis which asserts that neutrinos 
are massive particles, has been directly tested in exper- 



iments involving precise measurements of the endpoint 
of the beta decay spectrum, if one is interested in the 
mass eigenstates present in what is called electron neu- 
trino, or the kinematic behavior of the charged lepton 
produced in pion decay or tau decay, if one is interested 
in the mass eigenstates present in muon or tau neutri- 
nos, respectively. Nevertheless, such observations had 
generated so far only superior limits to the values of the 
neutrino masses and no absolute values of such quantities 
were measured. 

The second Pontecorvo's hypothesis, i.e., the mixing 
hypothesis, could be directly tested carefully analyzing 
the composition of a neutrino beam just after its creation, 
very close to its source. An ideal experiment would con- 
sist of positioning a detector sensitive to different neu- 
trino mass eigenstates very close to a neutrino source 
and observe its compositeness. Since detectors are only 
sensitive to neutrino flavor eingenstates, such an exper- 
iment cannot be realized. Nevertheless some hints on 
the neutrino compositeness could be achieved analyzing 
the flavor content of the recently created neutrino beam. 
The Pontecorvo mixing hypothesis foresee that close to 
their source neutrinos are found in a pure flavor eigen- 
state. Therefore, according to this Pontecorvo's hypothe- 
ses, very close to a reactor neutrinos have to be in pure 
electron flavor in the same way that very close to the 
pion decay pipe in an accelerator experiment, neutrinos 
have to be muon neutrinos or antineutrinos. 

Nevertheless, there are indications that this is not al- 
ways the case. Recent theoretical calculations of neutrino 
flux from nuclear reactors indicate that a larger than 
previously expected neutrino flux is produced |14L I15j . 
Such new fluxes are not entirely compatible with the 
short-baseline experiments which measure electron an- 
tineutrinos in distances of order 10 to 100 meters from 
nuclear reactors. Furthermore, the procedure of calibra- 
tion of the experiments GALLEX [16] and SAGE [17] 
which measured neutrinos within distances as small as 1 
m or so from the source raise also some incompatibility 
with observations and predicted neutrino flux according 
to these new theoretical calculations. Such incompati- 
bility has been called the anomaly of reactor antineu- 
trino [m [TS] and Gallium anomaly, respectively. Several 
different phenomena have been evoked to explain such 
anomalies [TH[T5] . 

In the present paper, we raise the possibility that the 
incompatibility of predictions and observations related to 
the reactor antineutrino and Gallium anomalies is a con- 
sequence of the usual interpretation of the Pontecorvo's 
quantum mixing hypothesis which defines, in a very fun- 
damental and unique way, what is the admixture of mass 
eigenstates in a specific flavor neutrino eigenstate. We 
will keep the usual interpretation that a neutrino pro- 
duced in a reaction in which a charged lepton is involved 
is a neutrino of the same flavor of this charged lepton. 
Therefore, the antineutrino produced in a /3-decay will 
be assumed to be of electron flavor. As well as, in a pion 
decay, once that a muon is involved, the corresponding 



neutrino will be assumed to be of muonic flavor and so 
on. Note that this is an arbitrary supposition once that 
neutrinos are not directed observed neither in creation 
nor in the detection processes. Nevertheless, different 
from the usual interpretation, we will assume that the 
admixture of mass eigenstates in the moment of neutrino 
creation is not unique but can vary for different neutri- 
nos produced in that reaction. This implies also that 
what is called an electron neutrino in the creation mo- 
ment can be a different combination of mass eigenstates 
from what is assumed to be an electron neutrino at the 
detection moment. Although unusual, we notice that 
such a hypothesis has never been tested so far and propi- 
tiates a possible explanation for the antineutrino reactor 
and Gallium anomalies, as we will see in the following. 
This new hypothesis and its consequences is what we call 
the Stochastic Neutrino Mixing Mechanism (SNMM). 

In order to appreciate how the SNMM works, we will 
analyze the particular case in which only two neutrinos 
are involved in the oscillation process. The generalization 
to the three neutrino case will be done in the next section. 
We propose relaxing the Pontecorvo's mixing hypothesis, 
allowing that neutrinos can be produced in an arbitrary 
superposition of neutrino mass eigenstates, each of them 
parametrized by a specific mixing angle 9c in the creation 
moment: 
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where 9c can assume, in principle, any value in the inter- 
val [0, ^]. The same assumption is made in the detection 
process, where the flavor state can also be identified in 
an arbitrary admixture of physical states, parametrized 
by a mixing angle at the detection moment in general 
different from the creation one, defined as 9d'- 
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Again, 9^ can assume any value in [0, f]. Under such as- 
sumption, after some distance L from the source to the 
detector, the Ve neutrino will present a survival probabil- 
ity calculated as 
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= cos^ {9d ''9c)- sin 29c sin 29 d sin^ ( — -^ ) , 

(3) 
where E is the neutrino energy and A?n^2 is the usual 
squared mass difference between the mass eigenstates in- 
volved in the oscillation process. Interesting enough, for 
the general case in which 9c ^ 9di this survival probabil- 
ity can be smaller than the unity even in short baselines 
in which L — > 0. Such behavior, which is not allowed in 
usual oscillation processes, is the essence of the solution 
of the Gallium and reactor neutrino anomalies which will 
be explored in the next section in the more realistic case 
envolving three neutrinos. 



II. THREE NEUTRINOS CASE AND THE 

SOLUTION TO THE GALLIUM AND REACTOR 

ANOMALIES 

We propose relaxing the Pontecorvo's mixing hypoth- 
esis, ahowing that each neutrino flavor eigenstate can be 
produced and detected in an arbitrary superposition of 
neutrino mass eigenstates around the usual admixture. 
In order to keep the success of neutrino oscillation ob- 
servations, we assume that neutrinos are created and 
detected most of the time around the usual superpo- 
sition of neutrino mass eigenstates which fit the oscil- 
lation phenomena, parametrized by the usual neutrino 
mixings sin^ 6*12 = 0.320 ± 0.050, sin^ 6*23 = 0.6131;°:°^^, 
and the recently measured sin^ 6'i3 = 0.025 ± 0.008, in 
3a [T3]. In general, these specific angles are going to be 
assumed only as the averaged values of the actual mix- 
ing angles. Under this simple assumption, we will con- 
clude that besides keeping the good fit of the observed 
long baseline neutrino oscillation phenomena, one can fit 
short baseline neutrino data setting a natural explana- 
tion for the anomaly of reactor antineutrino as well as 
Gallium anomalies. We assume, for simplicity, that such 
arbitrary superposition involves only the first two neu- 
trino families. Therefore only variations around 612 will 
be considered [TO] . 

The 3x3 mixing matrix at the moment of the neutrino 
creation (/7°) and at the detection moment (C/'') can be 
written as: 
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sin^y , d^'"- — cos 9c,d and s'^' — 
singed, and 0c,d can assume values in the interval [0, 7r/2]. 
The one particle electron neutrino survival probability 
can be computed: 
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(4) 
where 7 and /? run from 1 to 3. And then, averaging over 
different mixing angles, the total probability becomes: 
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where f{9c) and f{9d) are the distribution functions of 
the mixing angles involving only the electronic-muonic 
channel at the creation and detection instants, respec- 
tively. To keep the good fit of oscillation hypothesis with 
solar neutrino data and long baseline reactor observa- 
tions, we choose these distribution functions as: 
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which guarantees that mixing angles Oc4 will present an 
average value given by 612- In the above equation, ac,d 
are the Gaussian widths at the creation and detection 
instants, respectively, and we will assume, for simplicity, 
etc = Old = a. The normalization, A^c,d, is computed 

by imposing j^ f{dc,d)d9c,d = 1- Note that in the 
limit case when a -^ 0, we recover the usual Pontecorvo 
mechanism. 

Using all data of GALLEX and SAGE experiments [HI 
[T7] (see also Ref. gO|), old reactors [HI [H] as weU as 
the Daya Bay data jllj with a free normalization found 
according to the new flux calculations for reactor ex- 
periments, we perform a global analysis through the x^ 
method, defining: 

4 

where i and j correspond to each one of the four sets of 
experiments indicated by the labels appearing in Fig. IT] 
i,j ^ I for GALLEX and SAGE, i,j ^ 2 for old reac- 
tor experiments |21j . z,j =3 for Daya Bay and i,j = 4 
for Chooz and Palo Verde. Wij is the correlation matrix 
[2] , in which correlations between data coming from re- 
actors described by i,j = 2 are taken into account, while 
no correlation among other data is assumed, and col- 
umn vector i?^ collects the experimental data while i?* 
the corresponding theoretical predictions for reactor and 
gallium experiments. For reactors one has: 
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where S{E) is the energy neutrino spectrum which can be 
found in reference [22] and (t{E) is the cross section [14j . 
In GALLEX and SAGE radioactive calibration experi- 
ments, the reactions of electron capture produce neutri- 
nos of fixed energies. This implies: 
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gallium 



/ dVL-^ Y.^ {B.R.\cj,P,^^,^ {L, E,) 
JdVL-^j:^{B.R.)^a, 



(8) 



and the branching ratio (B.R.), the cross section (ctj) and 
the detector specifications are found in Tables 1 and 2 of 
Ref. [20] and references therein. 

The set of data includes 4 points from GALLEX and 
SAGE [5D], 21 from old reactors [2] as well as 6 from 
Daya Bay [TT]. We obtain the best fit value for a = 0.174 
varying in the intervals [0.141,0.201], [0.117,0.222] and 
[0.067,0.249] at 90, 95 and 99% C.L., respectively. 

This probability fits the data with a minimum 
Xmin — 39.08 which can be compared with the one 
obtained from the usual Pontecorvo's hypothesis result- 
ing x^ = 48.24, for 31 — 1 degrees of freedom. The best fit 
of SNMM as well as the fit coming from the usual Pon- 
tecorvo's hypothesis are shown in Fig. [T] where it can be 
seen that the SNMM provides a possible explanation for 
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FIG. 1: Comparison between the standard Pontecorvo's 
hypothesis prediction and the SNMM one plotted using 
the best fit value for the Gaussian width a = 0.174. 
The experimental points are distributed in the following 
way: 1. GALLEX and SAGE data; 2. old reactors [2T] : 
3. Daya Bay data with free normalization; 4. Palo Verde 
and CHOOZ. 
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FIG. 2: The Gaussian width a calculated to each 
set of experiments: GALLEX and SAGE, Reactors, 
LSND, Fermi and NuTeV. Points indicate the best fit 
at 68% C.L. and the curve shows the fitting to these 
points of the functional form a{E) = A+ {B/E)" , tak- 
ing the best fit parameters A = 0.012, B = 0.076 MeV 
and n = 0.565. 



short-baseline neutrino disappearance, something which 
is not allowed in usual neutrino oscillations based on Pon- 
tecorvo's original hypotheses. 

Before introducing our conclusions, we will add a possi- 
ble extension of the SNMM scenario. Up to now, we have 
discussed the relaxation of the Pontecorvo's hypotheses 
assuming a Gaussian distribution for the mixing angle 9i2 
characterized by a constant value of the corresponding 
Gaussian width a. Here we will observe that the SNMM 
can nicely fit several experiments assuming an energy 
dependence of such width. This is a consequence of the 
fact that, differently from low energy short-baseline ex- 
periments, high energy short-baseline experiments do not 
present any appearance or disappearance neutrino phe- 
nomenon. In fact, the neutrino disappearance is more 
intense in GALLEX and SAGE ^'^Ar and ^^Cr sources 
than in reactor experiments. In the first case the ratio 
R is smaller than unity nearly 14% j20jand in reactors 
R is lower than unity nearly 6% [Hj. Note however that 
the neutrino energy released in "^^Ar and ^''Cr sources 
have average value of 740 KeV while neutrino from re- 
actors possess a wide range of energy with a peak in 
3.6 MeV. Similar fact occurs in accelerator experiments. 
LSND [23] shows an excess of electronic neutrinos for en- 
ergies about 30 MeV. MiniBooNE |21] searched in two 
channels v^ — >■ v^ a-nd D^ — ^ i'e for oscillations. In the 
energy range of 200 < i?/MeV < 1250 was found sig- 
nal of oscillation in both channels, however, data sug- 
gest that the excess of events decreases when the neu- 
trino/antineutrino energies increase. The experiment de- 
scribed in Ref. |25 j. which we refer to as Fermi, worked in 
a different scale of energy, with peak in 30 GeV, searching 
for oscillation in the v^ — > Ve channel and did not report 
any signal of oscillations. The same happened in NuTeV 
experiment |26j . Executed with an average energy of 
about 200 GeV, it did not find signal of oscillations in 



both channels v,. 



and v.. 



The only possible 



exception is the experiment KARMEN [2^ that was exe- 
cuted with energies of about 15 MeV, lower than LSND. 
Although it did not find any compelling excess related 
to the background, its measurement was associated with 
large uncertainties. 

The above cited experiments suggest that there is a 
relation between appearance/disappearance phenomena 
with the energy. Identifying this possible dependence, 
we independently calculate the free parameter a for each 
one of the following groups of experiments: 1. GALLEX 
and SAGE [UlIT^, 2. all reactor data analyzed in [H] 
and Daya Bay [ll], 3. LSND [23], 4. Fermi [25] and 5. 
NuTeV [26], which result is indicated by the points and 
their uncertainties at 68% C.L. in Fig.[2| To fit all data, 
we propose that the width have an energy dependence 
a{E) = A+ (B/E)". Taking the best fit parameters 
{A = 0.012, B = 0.076 MeV and n ^ 0.565), we also 
show this curve in Fig. [2] 



III. CONCLUSIONS AND FINAL COMMENTS 

The SNMM can accommodate data that indicate dis- 
appearance of electronic neutrinos/antineutrinos in very 
short-baseline experiments. Assuming that neutrino mix- 
ing angles can vary according to a Gaussian distribu- 
tion around a preferable value given by the usual mix- 
ing angle ^12, a gaussian width a around 0.17 can fit 
all experimental data in what is called Gallium and an- 
tineutrino reactor anomalies. Furthermore, identifying 
an energy dependence in short-baseline neutrino disap- 
pearance/appearance phenomena, we could explain dif- 
ferent behaviors of high and low energy neutrino exper- 
iments assuming an energy dependence of the Gaussian 
width a which characterizes the SNMM. 

A few final comments are in order. First, we do not ex- 
pect significant modifications of previous analyzes involv- 



ing solar, atmospheric and other long-baseline neutrino 
experiments due to the implementation of the SNMM. 
Only modifications of few percents in the initial neutrino 
flux predictions as well in the detection rate calculations 
in the analyses of those experiments can appear due to 
the SNMM. They can be accommodated in several un- 
certainties present in these analyses and will not substan- 
tially alter their results [19j. 

Secondly, it is often assumed that contributions to neu- 
trino masses come from new physics, while neutrino in- 
teractions are given by the Standard Model. Neverthe- 
less, neutrinos are observed only indirectly through their 
interactions which produce charged leptons. This repre- 
sents a challenge to implement realistic neutrino sectors 
in any model describing this particle. Some previous ar- 
ticles propose new approaches. An interesting discus- 
sion about the definition of a fiavor neutrino state and 
its relation with physical neutrino states can be found 
in Ref. j28ri31) in which possible mechanisms which can 
generate non trivial mixing matrices that can be different 
in the neutrino creation and in the neutrino detection are 
discussed. This is one of the requirements to implement 
the SNMM and can inspire the proposition of models 
which accomplish the mechanism. 

We also propose possible tests to the SNMM. A niuon 
neutrino detector located near a reactor could exclude 
this hypothesis in case no muon neutrino would be found. 
Such a detector could be based on muon neutrino elastic 



scattering on electrons, in a similar way which is dis- 
cussed in Ref. [33], which is able to explore the muon 
neutrino at zero distance by v^e scattering. Also, ra- 
dioactive electron neutrino sources allocated inside ex- 
periments able to detect neutrinos through both charged 
and neutral currents channels (like as SNO |32j) would 
test the SNMM hypothesis. A non-oscillation effect in 
the neutral current measurement and an oscillation ef- 
fect in the charged current can favor SNMM in contrast 
to the sterile hypothesis, while a oscillation effect in the 
NC and CC measurement can indicate the presence of 
sterile neutrinos [15] . 

Finally, when we include an energy dependent gaussian 
width a we can fulfill the constraints on oscillation effects 
from low energy experiments, such as reactor and Gal- 
lium experiments, as well as the high energy experiments, 
such as FERMI, LSND, and NUTEV. We show in Fig. 2 
that a weak energy dependence is sufficient to achieve a 
nice and consistent picture of SNMM as a solution to the 
reactor and Gallium anomalies. 
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